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Symbols 

E  0  =  Vector  magnitude  of  incident  field 

B  0  =  Geodesic  field  vector 

)c0=  co/c  =  Free  space  wave  number 

H  '  =  Relative  permeability  constant 

c  '=  Complex  relative  dielectric  constant 

c  /=  Complex  relative  dielectric  tensor 

n  =  Complex  vector  index  of  refraction 

y-  J  k<>  n  =  Complex  vector  propagation  constant 
I         0=  Angle  between  y  and  B0 
I         £  =  Cos"1   ae  .  az 

X/  y,  Z  =  Cartesian  coordinates 

r,  d,(p=  Spherical  coordinates 

£=  Small  displacement  along  X-axis 

F  =  Force  vector 

e  =  Positive  magnitude  of  electron  charge 

m  =  Mass  of  electron 

-  e/m  =  -1.77  x  10U  coul/Kg 

v=  Average  number  of  electron-ion  collisions  per  second 

r  =  Radius  vector 

V  =  Vector  velocity  of  an  electron 

N  =  Density  of  electrons  per  cubic  meter 

jLi0  =  4  77  x  10       h/m  =  Permeability  of  free  space 

€0  «  J: x  10   '    fd/m  =  Permitivity  of  free  space 

3  6tt 
^  =./  Ne  2         =  Resonant  angular  frequency,  lossless  plasma 

03  =  ^B0    =  Cyclotron,  or  gyro,  frequency 
m 

E,  H  =  Electromagnetic  field  components 
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Symbols  (Continued) 

D  =  Electric  flux  density  vector 

t  =  Current  density  in  amps,  per  sq.  m. 

o  =  Surface  charge  density  in  coul.  per  sq.   m. 
■  s 

p  =  Charge  density  in  coul.  per  cubic  m. 
X  =  cop2/  CO2 

Y  =  o)b/co=  VyZ    +  YZ  +  Y2 

x      y       z 

Y  =  Y  cos  B,  Y    =  Y  sin  9  cos  <o,  Y  =  Y    =  Y  sin  ©  sin  <a=  Y  cos  ib 
x  y  z      L  r 

YT=V2+      2       =  V      2  _     2  =Ysin0 

x        y  L 

r  =  i-j  v/u 

M2  =  n2  -  1  =  n    .    n  -1 


a .  =  Unit  vectors 
i 

A.  =  a     x  a     =  a,  sin 
1        6        z         1 

C 

A0  =  (a    x    a     )     a. 
2          0         z          z 

=  a~2 

sin  £ 

X  =  1  -     r  X 

2         2 

r   -  y 

•N^ 

Y  =       X     Y 

r2  -x2 

Z=  1  -    X 

V 

• 

E        =  V  p2            2 
XY          E     +  E 

x      y 

xy 

*     -1  E 
=  tan          y 

E 
X 

0)    =  Critical  frequency 

CO      =  Ordinary  wave  critical  frequency 

co  T  =  Left  circularly  polarized  wave  critical  frequency 

co  D  -  Right  circularly  polarized  wave  critical  frequency 

CK 

<j  =  Lorentz  conductivity  tensor 
p  =  Lorentz  resistivity  tensor 
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Symbols  (Continued) 

I  =  Identity  tensor 

s  =  Length  of  path  through  the  ionosphere 

o 

c  p»  3  x  10     m/sec  =  Speed  of  light  in  a  vacuum 

V  =  Group  speed 

g 

V  =  Phase  speed 
P 

a.    =  Attenuation  constant 

]8    =  Phase  constant 

0O  =  Angle  of  incidence 

0    =  Angle  of  refraction 

R    =  Measure  of  polarization 

n..  =  n,  n.  =  Components  of  refractive  tensor 
ij         1    j 


1.     PLASMA  FREQUENCY 
1.1    Lossless  Case 

Y    ]jj//  ww.—  -   «~-~-  „..— . „  -  ™™,   «    ,m 

p„=""We4       of  free  electrons.     If  the  electrons  are  compressed 


D 


Consider  a  lossless  plasma  having  a  density  N  ( „§    ) 


■+ 


X       laterally  through  a  small  displacement,   £  ,  along  the 


x=0     v=  t 

X-axis,  a  surface  charge  density 

^^  -ps  =  -Ne  «  (1) 

may  be  considered  to  exist  at  y  =  £  .    Then,  from  the  definition  of  the 
electric  displacement  vector,  from  y=  o    to  x=  £  ' 

D    -p    .  (2) 

Hence,  D    =  *0  E    =  Ne  £    ,  (3) 

XX 


or  the  force  F     is  given  by,  _ 

X  F     =  _eE     =  -  Ne      '  £.  (4) 


X 

Since  force  is  given  by  the  product  of  mass  with  acceleration, 

2 
m     d     !    +    Ne2      4=0.  (5) 

dt2  C° 

2  2 

Substituting    to      =    Ne        into  (5),  the  harmonic  equation  for  free  oscill- 

m  c  0 

ations  of  the  electrons  is  obtained,  namely, 

14-  +  %  £=  0.  (6) 

dt*  p 

The  solution  of  (6)  is  sinusoidal  with  an  angular  frequency    co    .    Thus    co 

is  called  the  resonant  angular  frequency  for  the  lossless  plasma. 

1.  2    Lossy  Case 

Suppose  there  are  an  average  number,  v  ,  of  collisions  per  second 
between  electrons  and  positive  ions.    This  introduces  an  average  loss  of 
momentum,  or  an  equivalent  fluid  resistance  term,  into  equation  (5),  which 

equation  becomes, 

2  ? 

m     d     £       +  vm      d  £      +    Ne        t  =  q.  (7) 

dt2  dt  c0 

Pt 
Dividing  by  m,  and  using  the  operational  method  of  substituting  £=£0  e 


into  equation  (7),  the  determinantal  equation  is  obtained, 

(p2  +  VP+  to    2)  4=0.  (8) 

Equation  (8)  yields  the  characteristic  roots,  . 

P--I//2  ±jo)pV    i-  (1/(2  to^)z  (9) 

From  equation  (9),  it  may  be  concluded  that  ionic  collisions  in  the 
plasma  introduce  a  time  decay  constant  of  2/v,  and  reduce  the  plasma 
frequency  by  a  factor    V  j_  /  v/^        ,\2  .    For  a  low  loss  plasma,  the 
latter  factor  usually  is  ignored. 

2.     HOMOGENEOUS  ISOTROPIC  IONOSPHERE 
2.1    Low  Loss  Cas_e 

Consider  a  small  collision  frequency  for  a  homo- 


geneous isotropic  ionosphere  of  electron  density 

N(^L)   •     Let  a  uniform  plane  wave  front,  E0ge 
mo 

be  incident  upon  the  ionosphere  at  an  angle  of 


incidence  90.     If  r  is  a  radius  vector  describing  the 
path  of  motion  of  an  electron  having  a  velocity  V,  the  force  F  upon  an 
electron  is  given  by, 


or 


F  =  -eE  =  mdr       +ymdr 

at2            dt      ' 

(10) 

dV     +  vV    =  -  e     E. 

(11) 

dt  m 


Substituting  V  e  into  (11)  and  solving  for  the  steady  state  velocity, 

V  =  -  e    E  =      -  e  E 


,  .  (12) 

m  (jto+  V)         jm(  a>  -j  V) 

The  equivalent  current  density  due  to  the  motion  of  electrons  then  becomes, 

;=pV  =  -NeV    =Ne2E  _ 

m  {v+  ) 03) 

Equations  (12)  and  (13)  show  that  the  electron  path,  velocity,  and  current 

density  are  linear  in  the  direction  of  E. 

In  order  to  see  the  effects  upon  the  electromagnetic  wave,  Maxwell's 
equations  in  the  steady  state  phasor  form  will  be  used. 


The  Maxwell's  equations  are,  div  B  =  0,  div  D  =  0,  curl  E  =  -  _d b  / 

at 

curl  H  =  7  +  dp      ,  B  =  ju'  M  o  H  ,  5  =  e'  <r0  E  .  (14) 

St 
The  corresponding  complex  form  is, 


v  .  15=0,  v.E=0,  vxE  =  -jo)|Ll0H, 
V  x  H  =  i  +  jw  €0  E. 

Substituting  from  equation  (13), 

2  - 


(15) 


vx  H  =  Ne     E       +  jco  e0E  =  j<oe0  ^ 
m  (v  +  j  w) 

[1  +  - Ne2       ,       ]   E  , 


(16) 


or, 


jco  c0  m(y+  jco) 


VX  H  =    JO)  €0[  1  + 


cop2  (v-  jco) 
ju<i/z  +  CO2  )      _ 


E  = 


JCO  €, 


COr 


V  tu„ 


it  +  co2         coV  +  to2)    _ 


E  ,     (17) 


2  2 

in  which  to      =  Ne       was  also  substituted  in  the  right  member  of  equation 

m  c  0 
(17). 

From  equation  (17),  an  equivalent  complex  dielectric  constant  c      is 

obta  ined ,  «  0 


OJT 


<  =C  '-j  c"  =  (1- 

c  2,    2 

to  +  y 


J  -j 


_i     ^    CO 


co2+„2 


.       08) 


Returning  to  euqations  (15),  the  wave  equation  will  be  determined 

by  elimination. 

2  -  2  - 

Vx  (  vx  E)  sv(  V  .  E)  -  V     E  =  -V     E  =  -jcojU0* 

V  x  H  . 
Using  y   ■  &+  j  |3  =  jk0  V  c  '  ,  and  k0  =  oW  /!„  €0   =  co/c  ,  upon 

substituting  for  vx  H,  equation  (19)  becomes, 

V2  E    +  k2    e'  E  =  0. 

0      c 


(19) 


(20) 


For  y  ?»o,  n  =V  e'  will  be  evaluated. 

n-(c'-Jcl)  1/2='/" 


,  1/2 

(  H      *    ) 


2  c' 


-J-* 


2/T" 


(21) 


and  hence, 


y  =  j  k0n  M     k«  g//      +  j  k0jT~    . 
2/7T 


From  equation  (18)  with 

V* 

.0, 

y  fa 

vk 

a 

+ 

J  koV 

2A- 

to 

The  phase  velocity 

is 

V 
P 

=  to 

=  to/kn 

to 


(22) 


(23) 


^       V  1-  (  Wp)2  a/1-  (  Wp)2 


to 


to 


and  the  group  velocity  is , 


vg  =  _cUo_  =  c  71-  (wpj: 


d/3 


(24) 


(25) 


to 


2.2    Lossy  Case, 

For  the  lossy  case,  the  exact  value  of  n  =  JT~      should  be  used.    For 

^-""a  ?■        ^  —1 

this  purpose,  consider  /v/a  ±  .  b       ,        JS<^^\^'     '  B  =  tan      bA 


c^ 


sin  B/2  =4=-  Vl-jl 


77     VaTTP  . 


+  a 


7?  V"  a*  +  ba  -  a" 


Substituting  from  equation  (18)  into  (26), 

2 


n= 


to_ 


+  (1- 


2  +      2 
CO  1/ 


Ml-      V 


1/2 


1/2 


2  COp4  l^2 


2/     2^    2.2 

to   (to  +  v  ) 

2         to  4  u2 

)    +        P 


(26) 
1/2 

1/2 


O^2 


2/     2  j      2\2 
to   (  to    +  v    ) 


(27) 


2  j       2 
to    +  y 


Thus,  y  =  k0A$m  [n]and  j8  =  k0Re  [n].  (28) 

The  Re  [  n]  is  also  spoken  of  as  the  refractive  index  for  the  ordinary  wave. 

2.21    Critical  Frequency 

From  equation  (28),  the  phase  velocity  within  the  ionosphere  is, 

1 


V    = 


P      / 


/LioC'c  0        Re  [n] 


(29) 


Designating  the  angle  of  refraction  by  B    ,  from  Snell's  law, 

sin  8  =    sin  °° 

r       Re  [n]  (30) 

If  6r  <  ff/2,  the  wave  enters  the  ionosphere,  is  refracted,  and  passes 

through  the  ionosphere,  unless  the  right  member  of  (3  0)  is  greater  than 

unity. 


Figure  3 

If  the  right  member  of  (3  0)  is  greater  than  unity,  no  real  angle  0    exists  and 

the  wave  is  reflected.    Designating  the  critical  angular  frequency  by  a?      , 

this  frequency  is  determined  by  Q    =  tt/2.    Thus,  for  any  angle  of  incidence, 

r 


the  critical  frequency  is  determined  from, 

s  in  S0  =  Re  T  n  ]  . 
For  the  low  loss  case,  from  equations    (23)  and  (31), 

00  =  ^7 


sin 


05 


(  JJ2J' 


ttL 


or 


(31) 
(32) 
(33) 


o>c  =  cop    sec  90 

Equation  (33)  is  known  as  the  secant  law,  and  yields  05=0)     for 

c         p 

normal  incidence.    Substituting  the  numerical  values  of  e,  m,  and  c0  into 
(33), 

f   =  94Wsec  e0   /  (34) 

c 

in  which  f    is  in  cycles  per  second.     If  the  electron  density  is  given  per 

cubic  centimeter,  the  f    will  be  in  kilocycles  per  second. 

c 


At  a  given  angle  of  incidence,  frequencies  f  >  f  will  pass  through 
the  ionosphere,  and  frequencies  f  <  f      will  be  reflected. 

The  preceding  analysis  has  ignored  the  earth's  geodesic  field.  All 
waves  acting  accordingly  are  termed  ordinary  waves. 


3.     HOMOGENEOUS  ANISOTROPIC  IONOSPHERE 

3  .  0    Direction  Cosines 

In  order  to  take  into  consideration  the  effects  of  the  geodesic  field 

upon  an  electromagnetic  wave  propagating  in  the  ionosphere,  various 

vector  directions  must  be  taken  into  consideration.    Hence,  a  table  of 

direction  cosines  is  necessary.     Such  a  table  may  be  formulated  either 

by  projections  or  by  spherical  trigonometry,  or  by  a  combination  of 

both  methods.    Referring  to  the  figure: 

cos  j/)  =  cos  0  cost/2  +  sin  9  X 
sin  t/2  cos  (  tt/2  -<p) 

cos  0=  sin  9  sin  (p 

a    .  a     =  cos  j/f  sin  9  sin  tp 

X*;?    -  a  , .  a     =  sin  \b 

d^  0         Z  r 

-  a  .   .  a     =  -cos  lb  sin  tf 

lj)        y  , 

^    -  a  ,    .a    =  -cos  ilicos  9' 
0         X 

cos  0=  cos  j/jcos  tt/2  +  sin  j/j  x 
sin  tt/2  cos  9' 

cos  &'  =  cos  9 


sin  0 


sin  e'  =  Vl-.£os 


2  9 


sin     ij) 


Figure  4 


Vl-sin^Q  sin  2(p-cos2  eg  = 
'  sin2  j/) 

sin  9cos  (p 
sin  0 

sin  ij)=  Vl-  sin2  Qsin     tp 


From  the  above,  the  following  table  may  be  compiled: 


a 
X 

a 
y 

a 
z 

S 

ar 

cos  0 

sin  ©cos  <p 

s  in  6  s  in  cp 

0 

ae 

-sin  9 

COS    0COS    (p 

cos  9  sin  (p 

cos  Qsin  o 

sin  0 

a 

0 

-sin  <p 

COS  (p 

cos  (p 

sin  ^ 

0 

-sinScos    9s  in 

sin  ij) 

<P 

-sin^  gsincDCos  <p 

sin  0 

1 

sin  0 

(3  5) 


3  .1    Electron  Current  Density 

Assume  the  Z-axis  of  a  coordinate  system  to  lie  along  the  geodesic 
field  in  the  ionosphere,  and  assume  an  electromagnetic  field,  E  =  E0  aQ  , 
to  be  incident  at  angle  0O  upon  the  ionosphere.     Equation  (13)  gives  the 
current  density,  due  to  the  presence  of  the  electromagnetic  field  along 
with  the  free  electrons.     However,  it  ignores  the  presence  of  the  geodesic 
field,  the  geodesic  field  being  B0  a     . 

The  zero  order  component,   i0  of  equation  (13),  reacts  with  B    to 

produce  a  new  component  i.  ,  of  the  current  density,  normal  to  both 

i0  and  B    .     This  component    lies  within  the  XY  -  plane.     The  i    ,  component, 
z  1 

in  turn,  produces  another  component,   j_  ,  normal  to  both  l,  and  B     ,  and 

l  1  z 

hence  also  lies  within  the  XY  -  plane,  but  it  is  rotated  rf2  from  ti.    The 
new  component,    »      ,  produces  still  another  component,   L  /  normal  to 
I    within  the  XY  -  plane.     Since  \     has  been  rotated    it  from  f,   ,  it  is  in 

^  o  1 

the  same  direction  as  ^   /  but  of  the  opposite  sense.    Also,  the  component 
L^  ,  produced  by  t.  ,  and  is  also  of  the  opposite  sense.     Continuing,  an 
infinite  alternating  series  is  obtained  for  the  two  directions  lying  within 
the  XY  -  plane.    These  series  will  converge  under  certain  conditions. 
They  will  be  obtained  analytically  in  the  following    analysis. 

In  the  presence  of  the  geodesic  field,  the  force  equation  (11),  for  each 
electron,  must  be  modified  in  accordance  with  ampere's  force  law  to 
become, 


m 


d  V 
dt 


+  L»mV+eVxB    =  -eE  . 


(3  6) 


Equation  (3  6)  is  to  be  solved  by  iteration.    That  is,  the  solution  V0  of 
equation  (12)  will  be  substituted  into  e  V  x  B  and  a  correction  term  V.  , 
obtained.     The  V    term  will  then  be  substituted  into  e  V  x  B  to  obtain  another 
correction  term  V    ,  etc.     Before  substitution,  the  successive  cross  products 
will  be  formulated.     From  table  (35), 


Figure  5 


But  from  table  (3  5), 


a>a     sin  0+  a     cos  0cos  ©+  a     cos  0 

e     x  y  z 


sin  <p,  and 

A    =  a     xa     =  a.  sin  T  ,  T 
1  @        z         1  a      b 


(3  7) 


cos 


-1 


Al= 


Hence, 
and 


a  a  a 

x  y  z 

-  sin  0  cos  0cos  cp  cos  0  sin  <p 
0  0  1 


a    cos  0  cos  (0  +  a      sin  0 

x  y 


A0=  (  a^x  a    )  x  a     =  a0    sin  F  = 


(38) 


2      v     0         z'         z 

a  a  a 

X  7  z 

cos  0  cos  (psin  0  0 

0  0  1 


cos  0  cos  <p 


=  a    s  in  0  -a 
X  y 


(3  9) 


E     =  -E0  sin  8  ,  E     =  E.  cos  9 cos  tp,  E    = 

x  y  z 

E0  cos  0  sin  <p  . 


A,  =  a      E     -  a      E 

i      x    y      y     x 

A     =  -  (a      E    +  a      E   ) . 

2  x    x     y    y 


(40) 
(41) 
(42) 


From  equation  (38)  and  (39),  it  may  be  seen  that  both  A    and  A„    lie 
within  the  XY  -  plane,  and  that  their  slopes  are  negative  reciprocals.     That 

is  to  say,  A      is  rotated  tt/2    from  A..     Furthermore,  since  each  unit  vector 

a.   ,  j  =  1,  2  , ,  is  within  the  XY  -  plane  ,  each  a.       .  =  a .  x  a        also  is 

;.  i  +  1        i        z 


within  the  XY  -  plane,  and    a.       .  is  in  space  quadrature  with  a.. 

l+l  r  1 

Incidentally,  the  unit  vectors  a    ,  a     and  a fi  are  also  coplanar,  with 
a9  at  an  angle  t/2  +  £with  a    . 

For  the  iteration,  from  equations  (12)  and  (10), 

-    =  -  eE8  I  ft  p  =  e  voXb"=  jm  (W  -J||V..  (43) 

J  m  (  0)-  ji/)         1  z  * 

Thus,  upon  defining  the  gyro  frequency, 

ab  =  GBz  (44) 

m 

Vj=  e2  EoQ  x  Bz     =  -e  E0  o>b  Aj       = 


[jm  (  u)-]v)T      m  (  a)  -  ]v) 

-eg)     E0  sing     -  (45) 
/ i — 72"       * 


Repeating  the  procedure  for 


F     =-eV,xB    =  j  m  (o3-  ]v  )  V0    ,  (46) 

j  mz  (  oJ-'-jl/)3 

9 

-j  e  03,      E0    sin  £ 

£ a2 

m2  (  03  -  jvr 


For 


-, — TT 5    '     5         1 

(  03-  ]V)  b 


VA=~  J  e  «k  EoSinf    - 


7 
m  (  03-  jf) 


(47) 


F3  =  -e  V2x  Bz=  j  m  (  03-  jl>)  V      ,  (48) 

V3=  e  o>b3  E0sin  £ 

a.   ,  a0  =  ♦  a    .  (49) 

m  (  03-  jvr 

C°ntinUing'  v<-1.Mb4E..toC    -      S         "  (50) 

m  (03-  iv)b  4  '     4  Z 

V    =  -e  03,     E0  sin  f      -        -        -  /cl\ 

5  b      °  >       ar    ,  a     =a.      ,  (51) 


a6    'a6  =32   '  (52) 


The  above  procedure  may  be  continued  indefinitely. 


Upon  summing  the  corrective  components  (43)  through  (52)  ad  infinitum. 


factoring  -  j  to  e0,  and  substituting  to     ,  the  velocity  vector  for  an  electron 


Ne 
becomes, 

V  =  +  J  coe0E( 


co. 


Ne 

co  (  to-  j  i^6 
to  2    to  6 

_p b_ 

co(  to  -  )V) 


Q     to  (  to-  iv) 


+  j  a. 


co  2    w.  <*>   2  ^3 

p       ^b       +         P       a 


;in£   - 


sin  f , 


to2     to  2 


to(  to-  jf)z        to(  to-  jy) 
2  4 

k'r,  CO,    ^ 

+  P  b  + 


co  (  to-  )v)  3       to  (  to  -  jv)  5 


(53) 


or, 


V  =  +  j  toe  „E„    .  a 


to 


to  z     co 


Ne 


7 


+  j  a. 


(  co-  iv) 


co  (  to-  iv)' 


to. 


1  + 


(  co  -  iv)' 


CO    4 

+       b 


177x4 


(  co  -  iv) 

to  4 

b 

(  to-  ]v)' 


sinC  -a2       WP2    "b* 


sin 


1. 


w(  a)  -  ii/)3 


1+      "b 


+ 


(  co  -  j  v) 


(54) 


.J 


The  infinite  series  within  the  brackets  of  equation  (54)  converges  for 

<  1    ,  (55) 


cob2 


2  2 

to    +  v 


that  is,  for 


co> 


to,     VI- 


(  •  cob)' 


Hence,  for  the  frequency  range  of  convergence,  the  brackets  become, 


(  co-  j  v)2 


=  (  co-  )v)' 


(56) 


(57) 


1  -  co  2/  (  to-  iv)       (co-  iv)    -  CO 


2     ' 


and  hence, 


V=    +  j  cognEn 
Ne 


to 


tOu  co      sin  f 

+  b    WP  S  y 

~1TZ  ...  2  - 


to(co-jv)      co[(  co-  ]vr  -  to^] 


jara2 


CO 


co  -  iv 


(58) 
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3.2    Path  of  An  Electron 

For  an  examination  of  the  path  followed  by  an  electron,  it  is  preferable 
to  eliminate  to  2  in  equation  (58),  and  integrate  in  time  to  obtain  a  radius 
vector,  r  ,  whose  terminus  describes  the  path.    Accordingly, 


1 +     ^b/l-  cos_2  8sin  2 


r    =    E0     -z-l  an       _J +     wb,/  1-  cos_    9  sin      <p 

0       m 


u(  a)-  jv)  2  2  _ 

CO    [  (   C0-  J    V)       -    03)3       J 


]«,-a 


"b 


to-  j  y 
in  which,  from  equation  (37),  sin  >*     has  been  substituted  in  the  form 


(59) 


r~       2      2 

sin  £  =   v   1  -  cos       8  sin    p     .  (60) 

Assuming  c«  0,  equation  (59)  shows  that  as  the  electron  attempts  to 
vibrate  in  a  path  parallel  with  the  exciting  field,  it  has  an  additional 
elliptical  component  within  a  plane  normal  to  the  geodesic  field,  and 
resonates  at  to  =  to    .    For  frequencies  to<  to^  /  the  ionosphere  acts  as 
a  conducting  medium, for  the  current  density  produced  by  the  magnetic 
field,  rather  than  as  a  dielectric  medium. 

The  electron  path,  as  it  spirals  about  the  direction  of  the  exciting 
electric  field,  is  somewhat  like  a  trochoidal  epicycloid.    The  electron 
current  density,   i  =  -  NeV,  may  be  considered  as  a  source  of  radiation, 
a  vector  potential  formulated,  and  the  resulting  fields  computed  there- 
from in  accordance  with  Huygen's  principle.     Thus,  the  electric  fields 
are  no  longer  confined  to  the  path  of  the  driving  field.     It  will  be  shown 
later  that  D,  H,  and  y  constitute  an  orthogonal  system,  with  E  not 
necessarily  orthogonal  to  y  . 

3.3    The  Curl  H  Equation 

The  current  density  will  be  formulated  from  equation  (58),  by  mul- 
tiplication with  -  Ne,  and  substituted  into  curl  H  in  order  to  formulate 
the  dielectric  tensor.  The  dielectric  tensor  subsequently  will  be  used 
to  determine  the  indices  of  refraction.    The  resulting  curl  equation    is, 
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V  x  H  =  j  wc0EJa, 


1-  wp 


a)  (oj  -  ji>) 


-  j  a. 


0),     0)  „     sin  £ 
2  2 

03  [(03~  )V  )      "   0)5    ] 


5, 


^b    V  sin  g 


2  2 

0)(0)-  j^)[  (  03-  jl^)      "  o^    ] 


(61) 


The  wave  corresponding  to  the  first  term  of  the  right  member  of  equation 
(61)  is  sometimes  called  the  ordinary  wave,  and  the  wave  corresponding  to 
the  second  and  third  terms  is  companionably  called  the  extraordinary  wave. 


3.31    Cartesian  Coordinates 

Using  E    =  a     E     +a     E      +a     E    along  with  equations  (41)  and  (42), 

0      x    x      y    y       zz 

equation  (61)  may  be  re-written, 


i v  x  H  =  a 

j    03Co  XI      X 


-E      ^l 

X 


1  - 


03. 


O)   (  03-  )V) 


-JE        ^b    "p 

y  — 


0) 


(03-  )V)      -0^ 


9  2 

U>(  0)  -  )V)  [  (  03-  \v)    -  o>b  ] 


v  a    JE 

y-f  y 


03 
1-  P 


03(0J  ~  )V  ) 


03,       03    2 

J  E  b       p 

X  ' 


-  E. 


03,  2      0) 

b         p 


2  2 

O3[(or  ji^    -  wb   ] 


2  o 

u>(u)  -  ji>)  [(ov- ji>)    -  03b   ] 


+   a    E 
z    z 


ii-  V     1 

03    (03-  )V)J 


Now  let 


X   =  1  - 


03. 


03,  03 

b         P 


2. 


2  2 

0)(03-  JV)  0)(03  -  J  I/)    [  (  03-  jv)      -  03b      ] 


1  - 


0)    2    (03  -  jv) 

" 9 1- 

0)   [  (  03-  JV)   -«qJB 

2 


Y<=         "b     "p 


0)[(  03-  )V)2   -   03b2] 


,Z     =1- 


03 


0)   (03  -  jy  ) 


and  substitute  into  (62)  to  obtain, 

—        V  x  H  =  a     (E    X*jE     ?)  +  a     (E     X+jE     Y)  +  a     E     Z 

j  o)c0  xx  y  yy  x  zz 


(62) 


(63) 
(64) 


(65) 
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For  cross  reference  to  literature  employing  the  standard  URSI  symbols, 
X=oj2  Y=^h         r  =  1  "  J    v/  0)    i  (66) 

0)Z  to 

the  symbols  in  (63)  and  (64)  are  equivalent  to, 

X  =  1-  X  r  ,       Y  =     X    Y  ,  Z  =  1  -       X        .  (67) 

2         2  2         2 

r   -  y  r   -  y  r 

The  symbols  (67)  will  be  used  later,  but  for  the  present  it  is  more  conven- 
ient to  retain  the  symbols  (63)  and  (64). 

3  .32    Rotating  Coordinates 

The  a      and  a      components  of  equation  (65)  will  be  broken  into  left 
X  Y 

and  right  hand  circular  components.    Re-arranging  equation  (65), 

1 V  x  H  =    a     (E     X-jE    Y  )  +  j  a     (EY-jE    X )  +  a     E     Z  . 

"7^77  xx  y  y    x  y         z   z      (68) 

or, 

J V  x  H  =    E    (  a      X  +  j  a     Y  )  -  j  E     (a     Y+ja     X).  (69) 

jcoc0  xx  y  y     x  y 

Now  let, 

a     X  +  ja    Y'=A.  (a     +ja    )  +  A0  (  a    -ja)     ,  (70) 

x  y        i     x        y       2      x       y 

and  solve  for  the  undetermined  constants  A.  and  A     , 


~  Y     A-     V 

A,  +  A„  =  X  A    = 


1        2  1  2 

A    -A     =Y  A     =  I    "    Y 

1  I  L  9 


(71) 


also  let, 


ax       jay5r=Bi(5x  +  iay)  +  B2(ax"iay  '  '  (72) 


and  obtain, 


(73) 


Bl-B2=X  B2=X^  =  -A2 

Substituting  from  (71)  and  (73)  into  equation  (69), 
1 V  x  H  =  (E    -  j  E   )  *    +    Y   (a     +  j  a    )+(E    +jE)x 

j  w€0  x       y      2         x        y        x       y 

X    -  Y   (a      -  j  a    )  +  (l     E     Z.  (74) 

T~       x        y        z    z 

z  13 


From  equations  (63)  and  (64),  2 

Y  =  i-    ty  (M-  M*  (Ufa  ^P2     =  1  -  ^p     ^~  J">  T   ^b3 
C0[  (  co-  jy)     -  CO  2] 


X  ±Y 


1  -       X 


r  ±  Y 
upon  letting , 
E 


co[(co-  ji>)  +  o)b][(o>-  jv)  -  ovl 

(75) 


2  2  -IE 

1/2  VE    +  E  ,0         =  tan         _v_ 

XV  X        Y  XV  E 


(76) 


equation  (74)  may  be  written, 


V  x  H  =  j  co  € 


r 


E        e 

XV 


-  J  0 


xy 


L 


(1-   J< — )  (a     +  ja)  +  E 

r  +  y      x       y 


XV 


XV 


(1  -      X     )  (a     -  j  a))  +   a      E     (  1  -   X_  ) 

T^T      x        v        *    *         r 


(77) 


3.321    The  Complex  Dielectric  Tensor 

Designating  the  left  hand  rotating  unit  vector  by  a     and  the  right  hand 

vector  by   a     ,  that  is  , 
R 


aT  =  a    +  j  a 

L     x        y 


a_   =  a     -  j  a 

R      x       y 


(78) 


and  letting , 


e/=l-_x ,    e'=l-_x ,    g>  l  -     x 

r+Y  r-Y  r 


equation  (77)  becomes, 

V  X  H  =  j   CO  €0 


0       %      0 


0        0 


-   j    0 

E         e         Xy 

xy 

E  ej0Xy 

xy 


(79) 


~        -1 

aL 

• 

*R 

a 
z 

U»           s—l 

=  JC0D.       (80) 


3.322    The  Complex  Indices  of  Refraction 

Equation  (80)  yields  three  characteristic  waves  propagating  in  a 
homogeneous  anisotropic  ionosphere.     There  are  two  circularly  polarized 
waves  within  the  XY  -  plane  rotating  in  opposite  senses,  and  a  linear 
wave  parallel  to  the  Z-axis.    These  waves  have  three  distinct  indices 
of  refraction. 
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Rationalizing  the  denominators  of  equation  (79), 


=  r{.  %  (<*-  wb) 


03r 


60 


2+  ^ 


-  J 


03 

1/ p 


03(o>^  +   l^)_ 


2         2 
0)[(03  +  cob)     +  ^    J 

V  I*  -  "J 

2         2 

03   [(o>  "   03b)      +  1/     ] 


K     ^P 


0>[(o>  +   U5b)      +   y2] 

2 
y     ^p 

2         2 

w  [  (oj  -  ov)    +  v  J 


(81) 
(82) 

(83) 


The  corresponding  indices  of  refraction  will  be  obtained  by  extracting  the 
square  roots  of  equations  (81),  (82)  and  (83).    The  ordinary  index  of  refrac- 
tion, n0  ,  is  identical  with  that  of  equation  (27).     Upon  applying  equations 
(26)  to  equations  (81)  and  (82),  the  left  and  right  hand  indices  become, 


n 


+ 


R  IL      to[(a>  ±  cobr  +  v    ]  -J      VL      W[(o)  ±  w 

"1  2  \  1/2 


wb)' 


J2  +  »2]J 


2       2 
CO  l^ 

P 


w[  (to  io^)2  +  i^] 
2 


1/2 


03    2    (03  ±    03b) 


03  r(o3  ±  03b)    +   V  2] 


/Ti-  Vto*«b2        j 2  +  J     V  " 

\(_       03  [(OJ  ±  03b)2  +   y     ]j  j      U)[(05±COb)      +  V 


1/2 
(84) 


It  should  be  kept  in  mind  that  equation  (80)  is  normalized,  and  hence 
indices  (84)  apply  only  for  propagation  of  these  characteristic  waves.    For 
wave  combinations  other  than  the  characteristic  waves,  there  will  be 
coupling  elements  within  the  dielectric  tensor,  and  a  dispersion  equation 
will  be  required  for  determining  the  indices.    Which  of  these  waves  appear 
depends  upon  the  direction  of  incidence  as  well  as  upon  the  frequency. 

3.323    Physical  Interpretation 

The  critical  angles  of  reflection,  corresponding  to  the  real  parts  of 
the  indices  of  refraction  in  equation  (84),  may  be  determined  by  equation 
(3  0).    This  will  be  done  for  a  low  loss  ionosphere,   v  &  0,  and  a  physical 
interpretation  formulated. 
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Setting  v  =  0  in  (84)  and  squaring,  is  equivalent  to  setting  v  =  0  in 
equations  (82)  and  (83).     Hence,  for  applying  equation  (30),  set 
2  2 


1  - 2_ 


(Ji  (  tO  ±    40b) 


=  sin      90  , 


°r'    a)     ±   0>b   W  -  to  2    sec     a0=  °    • 


(85) 


(86) 


Taking  the  positive  sign  in  equation  (86)  for  to  T   , 
to      =   -     ^   +  V  (  a)n  sec  e  0)2    +  (  to 


cL 


W  • 


(87) 


and  taking  the  negative  sign  for  to 


cR  ' 


«_„    =  +   Ji_  +  V  (to„sec  ©;)2  +  (  t*b)2      . 


'cR 


\7 


(88) 


Since  to       determined  by  equation  (81)  is  identical  with  to     determined 
by  equation  (33), 

wc     =  ^c    =    Wp  S6C    6o  ' 


while , 


/      2 

WcL  =  Va5c    + 


w 


(Ji 


b    , 


and 


WcR  =  *  Wc2  +  K^    +  _% 


2  2 

For  frequencies  sufficiently  high  such  that  (tOc)     >>  ((Ob*' 


Hence, 
(0 


40     x2 


I2  H 


r  «  to    +         b 
cL  c 


1/2  =  - 


wb 


:tf 


i+  »h. 


8«c2 


W    tO 


8   to  2 

c 


(89) 


=  to    +^b_+  —     .    (90) 
c    8  to 


(91) 


and, 


to  „  m  to    +  _^b     +    ^b      »  to      +    ^b 

cR         c      — — c         —7- 

8  tO           ^  2 
c 


(92) 


Since  the  three  characteristic  waves  have  distinct  indices  of  refraction, 
they  propagate  within  the  ionosphere  along  different  ray  paths,  and  have 
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different  attenuations.    The  attenuations  spoken  of  above  are  those  due 
to  $m  [n]  .    There  are  other  attenuations  due  to  numerous  anamolies  of 
the  ionosphere.     Because  of  the  different  critical  frequencies,  the  right 
hand  wave  may  be  reflected  at  a  frequency  for  which  the  other  waves 
pass  through  the  ionosphere. 

Recalling  that  the  series  in  equation  (54)  diverges  for  frequencies 
less  than  the  cyclotron  frequency  (or  gyro  frequency),  it  may  be  con- 
cluded that  for  a)  <  W, 


the  extraordinary  wave  confronts  a  conducting  boundary  and  is  reflected 

as  in  path  c    (Fig.   6).     For  frequencies  such  that  o),    <   05  <  05      /  the 

-~  be 

wave  enters  the  ionosphere,  suffers  a  refraction,  and  is  refracted  back 

to  the  earth,  as  in  path  b  (Fig.  6).     On  the  other  hand,  for  frequencies 

OJ  >    CO  D  /  the  wave  is  slightly  refracted  but  passes  through  the 
cR 

ionosphere  as  in  path  a_    (Fig.  6).    The  three  waves  do  not  necessarily 
all  exist  simultaneously,  depending  upon  the  angle  of  the  ray  path 
with  the  geodesic  field. 

3  .4    Azimuth  Component  of  Exciting  Field 

In  the  preceding  analysis,  the  exciting  field  was  considered 
vertically  polarized.    However,  a  wave  incident  upon  the  ionosphere 
from  a  horizontally  polarized  antenna  will  have  a  horizontal  component 
in  addition  to  a  vertically  polarized  component.     Hence,  it  becomes 
necessary  to  consider  the  azimuth  component  of  an  exciting  field. 
Re-writing  equation  (3  6)  for  convenience, 
F  =    -e  E  =  m  d  V    +  ymV+eVxB    ,     7    =    -NeV  , 

and  considering  the  exciting  field  in  the  form, 


(36) 
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E     =  En    (- a      sin  (p  +  a      cos  m)  , 
<P  X  z 

the  iteration  follows  analogously  to  that  for  E    .    Accordingly, 

8 

following  steps  (43),  etc., 

V0  =  -  E0  e    a  <p       ^  -     =  _  ey     x  B  =  j  m  (co  -  j^)  V 


j  m  (or  )v) 


1 


1    ' 


V,  =  e2  Eng(PxBz  =  e  E0  a5b  (-a    sin  <p+  az  cos  o)  x  az  = 

[j  m(to-  jy)]  m  (co  -  jy)z 

e  E0  co,    sin  <p  a 

a —  , 

m  (to  -  jv)2 


-  e     E„  B     w,    sinaa     x  a 


V 


j    mz  (  a?  -  j  vY 
,2  p   r     ,.,  2 


V     =      j  e     E°BZ  ^b     Sin  ^ 

0 


V4  = 


j  m2  (  co  -  j^  )' 


^ 


a    x  a 


■jeE  0  tOb     sin  <pai 


m 


(  to  -  jv)3 


eE0  co,3  sin  <p  a 


m  (  to  -  jv)4 
-e2  E0  B^  coh3  sin  cpa^xa     _     -j  e  E0  cob     sin  <pa 


X   » 


j    m  z  (  co  -  jy)5 


m  (co  -  jv)b 


X. 


(93) 


(94) 


(95) 

(96) 
(97) 

(98) 


Substituting  the  above  into  curl  H  ,  using  i  =  -NeV    ,   and 
considering  E0  a    =  a__  E._  +  a_  E_  ,  curl  H  = 


<P 


y    y 


j  W  ej   E0a^[l- 


.,  2 


z      z 


co  (to  -  jl>) 


]  -  ja    E 
X   V 


<V   F   "b  + 

oi     [(to  -  ii/)2 


to. 


U-\ 


+  -—      -  a    E 

(co  -  Jv)4       (to  -  W  )6 


+     % 


to, 


to  (to  -  jy)3 


(99) 


(to  -  jy)5      (to  -  jv)y 
Summing  the  series  over  the  range  of  frequencies  for  which  the  series 
converges , 


curl   H  =  j  40€0<{    a    E0 


1 


(oj  -  jv)2 


2  2 

to  .-    -      CO          to, 

l-_£ -jayEy-? ^ 

to(  to-  ]v)  _|  A        to(to-  jy 


-  a     E 

y    y 


p 


lh  (m  -  ]vY 


1-     wb 


b 

-  ivF 

l  "» 

"2 


(to-  j^)Z 


(100) 
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or, 


1  a     E0  [1-^E 

w  (a?-  jv) 


]   -jE      aP      ab 


C  a.,  -  j  aT 


03  -  jv) 


]  V. 


fioi) 


Substituting  for  a     E0  of  equation  (93),  and  using  E    and  E     , 


curl  H  =  j  o?€, 


-J  a 


E      03        o\ 

V       P         b 


03. 


03[(  05  ~  ]V)      ~    0?b     ] 


+  a     E    J  1-    -p 

y  y7    ;tut^7) 


2  o 


+  a     E      f  1- 
z     z    u 


05 


03  (03  -   jy)' 


(102) 


Let 


Substituting  from  equations  (63)  and  (64), 

curl  H=jo)€0/-jE       (a     Y+ja     X)  +  a      E     Z 

]  y       x  y  z     z 

a      Y+ja     X  =  A.  (  a     +ja    )  +  A„  (  a      -ja     ), 

x  y         i      x        y        2      x        y 


(103) 
(104) 


A    =  1/2  (Y    +  X) 
A2  -  1/2  ff    -  X)  . 


and  determine  the  arbitrary  coefficients  A.  and  A_   , 
A1  +  A2=Y 

Aj  -  A2  =  X      , 

Thus,   (104)  becomes, 

- —       curl  H  =  -  j  E     (X    +    Y)(a+ja)+jE      (X    -    ?)  / 

j  w€0  _y_  x       y        _l 

2  2 

(a    -  j  a    )    +  a     E      Z      . 
X         y  z    z 


(105) 


(106) 


Changing  to  the  URSI  symbols  by  substituting  from  equations  (67),   (75) 
and  (78), 


curl   H=  j  03cJ  -  j    E      (1-     X     )    a_   +  j    E    (1-     X       )  a_   + 

-f        r  + Y      L         -f      T  -Y      R 

E     (1-    X_)    a 

z       r       z 

Using  equations  (79) ,_ 


(107) 


Curl  H  =  j  o)  €( 


/ 

0 

0 

-j    E/2 

~h~ 

0 

€R 

0 

• 

J    E/2 

• 

SR 

0 

0 

1 

Ez 

a 
z_ 

=  j   03  D      . 


a  08) 
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Hence,  the  same  characteristic  waves  appear  as  in  the  case  of  the 
vertically  polarized  exciting    wave.    Which  of  these  waves  appear 
depends  upon  the  angle  of  incidence  and  direction  of  incidence, 
as  well  as  upon  the  frequency. 

3  .5    Horizontally  polarized  exciting  field 

Suppose  an  exciting  electric  field  is  incident  upon  the  ionosphere 
from  a  horizontally  polarized  antenna,  the  incident  field  being  in  the 
form , 

p  -         v        a    j  w  t  "  y  '  r  a 

L  ~    "±o_    e  a0      ,  (109) 

r 

-  a  ,    being  shown  in  figure  4.     It  may  be  seen  that  -  a  ,    is  composed  of 
two  orthogonal  components  in  the  a     and  a     directions. 

If,  in  equations  (40),  E0  is  replaced  with  Eofl  ,  and  if  in  equation 
(93),  EQ  is  replaced  with  E.      ,  then  from  table  (35),  EQof  equation  (109) 
may  be  expressed  by, 

0       8      8        <p      <p 


<-E0  a    =a_  E0fl  +  a     E.  (110) 


WlthE      =E      COS  9  Sin  <°  cos  Q 

°e  V  1-  sin2  9  sinz  <p     '  %"  E°  7T-  sin2  0  sin2  <p 

Also,  inequations  (40), 

/  p2   .  F  2    =  Eo  r         2^.2 

v    hv,    +  L„  v  1-  cos      9  sin 


(111) 
cos  "  8  sin'2'  <p     ,  (112) 


and, 

tan  0         =     y       =    -  cot  9    cos  <p  .  (113) 

XY      E 

X 

Therefore,  for  the  horizontally  polarized  exciting  electric  field  vector, 

equations  (77)  and  (107)  may  be  combined  intd  ^___ ' 

curl  H  =  jwc0     -Tl/2  [  e,]  \.Y~  f£   J^  e       XY  eR   aR  ]  J  I-  sin29cos2  93 
+  a    €q  cos  8  sin  <p    I     cos  8     sin  <p    E0 

!  yi-sin28sin2cp  <114> 

+  juxA  j  1/2  [cL'a*L  -  eR   aR  ]  sin  <p+  a"z   c0'cos  <pl    f  OS  *  ^ 

L  Jyi-  sin     8sin^<p 

To  determine  the  critical  frequencies  involved  in  equation  (114),  it  would 
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seem  preferable  to  consider  individually  the  various  normalized  components. 

4.    THE  DISPERSION  EQUATION 

4.1    Maxwell's  Equations  in  Complex  Form 

The  instantaneous  Maxwell  equations, 

div  B  =  0,  div  5=  p  ,  curl  E    =  -  a  B    ,  curl  H  =  d  5  m_. 

at  at      '  {     ' 

t    =a  E    =pV     =-  NeV  ,  5  =  €  '  €  0  E    ,  B    =  jli  '  jli  0  H    , 

are  possibly  more  useful  for  steady  state  time  harmonic  cases  when  con- 
sidered in  the  form, 

div  B  =  Ojdiv  5  =  0,  curl  E  =  -  j  OJ  /LlH,  curl  H  =  i  +  J  05  c  E  .  (116) 

To  show  that  div  E)  =  0,  first  formulate  the  equation  for  the  continuity 
of  charge  by  taking  div  curl  H  , 

div  curl  H    -  0  =  div  i  +  d_  div  D  =  div  i  +  d  p 


or 


at  at     ■  (117) 

div  t    =  -  an 

IT .  (118> 

But,  for  conductors, 

div  i  =  crdiv  E  =  cr      div  D  =jr  p    ,  (119) 

±£_+2_   p=    0  (120) 

at       € 

The  solution  of  (12  0)  may  be  written 

p=p0e'lf/^t      .  (121) 

) 
From  equation  (121),  a  charge  density  within  a  lossy  medium  must  vanish 

with  time,  and  hence  must  vanish  within  a  steady  state  condition.    Of 

course,  since  a  dielectric  is  considered  to  be  free  of  charges,  div  D  also 

vanishes  within  a  dielectric. 

Now  consider  the  case  of  a  spherical    wave  front  given  by, 

:;(fl      ,  j  0)t-y   .  r  (122) 

E  -  E0(8/  <p)      e, 

r 
which  is  sufficiently  remote  from  the  source  for, 


a  E 


0 


Rrf 


o  ,   a  e0  «  o  , 


ae  a<p 
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to  hold  locally.    Also  assume  y  =  a    y  .     Then, 

div  (  E0     e"y  '  *  )  =       1  (j_  ^ 

r  rz  sin  0/  Br 


EiLr_sin_8  e- y  .  r     j 


=  a  .    E, 


i_VE0e 


-  yr  e 
-y.  r 


-y.r  -y.r/  -     p  -y.r 

'  +er  £    =  -  y  .  Lo     ~      x  + 


y  •  H°_e 
r 


(123) 


But,  the  wave  is  assumed  to  be  a  radiated  wave  such  that  the  inverse 
square  terms  vanish.     Hence,  in  practice,  for  a  radiated  wave, 


div 


E0(6,cp)        -y  -r    =    -y.     e„  (  9  .  <P)      e    "V    •  r 


Likewise,  consider 

curl  (J,     e  "y  *  r)  =  _4_ 


r     sin  9 


r 

ydr 

E^J-y.  r 
r     e 


r  a 
6 

r,Eo9    e"  y  ■  r" 


(124) 


cos  0 


T     : — T    E0cp  e 
=  a    rc  sin  9 
r 


-y.  r 


r  sin  9    a 

a /dp 


<p 


r  sin  9  r       e 

r  ^ 


-y.r 


-a       r  sin  9 
9  E 

r     sin  0         ^ 


(  -ye  "y  •  r  )  + 


<p 


£f9_    (  -ye 


)      . 


(125) 


Again,  assuming  a  radiated  field  for  which  the  inverse  square  terms^vanish, 

curl  (E^_e  ~y  'T^  - 
r 


-y.r). 


)  =  ar  0  -  ag    (  -  y   E^     ~7    •  r)    +  g    (  -  y  E,  9 

r  r 


(126) 


But, 


yxE0  = 


y 

E„ 


<P 


°9      E° 


<P 


=  a    0-  a^yE.     +  a     y    E0rt 
r  0°<0       <p  °6 


(127) 
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Therefore,  in  practice, 

-y  .  r 


curl 


pL  (6  •  a 


=  -y    x    E0   (9    ,  <p) 

r 


-y  .  r 


(128) 
(129) 


Also,  in  case  the  relative  constant  is  a  tensor,  c',  then 

5  =  €0F'  .  E  . 

Therefore,  the  complex  form  of  Maxwell's  equations  for  radiated  fields 
may  be  written, 

y  .  H  =  0,  y  .   €o€'.  E  =  ,0,  7  =  -   NeV,  yxf-j  w^0H,    -  y  x  H  = 

-NeV  +  j  toe0  E  =  j  a)  €07'    .  E  ,  (130) 

in  which,  c  'is  defined  by 

"e  '  .  E    =  -  NeV/  joK0  +  E     ,  (131) 


The  wave  equation  becomes, 

y"x(yxE)-k20  ?.     2  =  0. 


4.11    Orientation  of  field  Vectors 


(132) 


From  equations  (129)  and  (130),  the  relative  orientation  of  the  field 
vectors  may  be  obtained.    For  this  purpose,  consider 

y    .  D    =  0,"y  .     H  =  0    ,  (133) 

from  which  it  may  be  concluded    that  y  is  normal  to  both  D  and  H. 
Then  consider, 

-yxH  =  +ju>5,  (134) 

from  which  it  appears  that  D  is  normal  to  both  y  and  H.    Therefore,  the 
three  vectors  D,  H,  and  y  constitute  an  orthogonal  triple.    Finally 
consider 

yxE  =  jo)(j0    H     ,  (135) 

from  which  it  is  seen  that  H  is  also  normal  to  the  plane  containing  yand  E. 

Assuming  y    to  be^t  an  arbitrary  angle  $  with  the  geodesic  field  B0  , 
the  relative  orientation  of  the  unit  vectors  is  then  illustrated  in  figure  7. 


Figure  7 
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Thus,  the  TEM  incident  wave  becomes  a  complex  TM  wave  within  the 
ionosphere. 


4.2    The  Dispersion  Equation 

For  y    =  j  k0  n  ,  the  wave  equation  reduces  to, 

nx(nxE)+e'.E=0.  (136) 

upon  expanding  the  cross  products  and  providing  for  all  components  of  the 
dielectric  tensor,  equation  (133)  may  be  expressed  as, 


/     2  2. 

-(n       +  n     ) 
y  z 


n      n 
y      X 


n      n 
z       X 


n      n 

x    y 


n       n 
X      z 


-(n2  +  n2)       n      n 
X        z  y      z 


n      n 
z      y 


i  2  ^       2. 
-(n    +  n     ) 

x      y 


'e"1 

X 


y 


"e  1 

X 

• 

E     4 

y  T 

E 
z 

y      w 

XX 

xy 

/ 

€yx 

yy 

€ 

z  X 

zy 

XX 


yz 


zz 


=  0 


(137) 


Combining  the  tensors, 


€  '       -n2_n2^       +nn 

xx       y       z      xy      x   y 


<^ 


+n     n 
y    X 


y 

c'         +n     n 
z  X         z      X 


e'     -n2  -n2 

yy       x     z 


e      +  n    n 
zy        z     y 


c        +  n    n 
Xz         X    z 

c'      +  n     n 
yz         y     z 

tf/      -n   2  -n 


zz 


y 


=  o 


(138) 


Equations  (136)  constitute  three  homogeneous  linear  equations  defin- 
ing the  components  of  E.  A  necessary  and  sufficient  condition  that  a 
solution  other  than  the  trivial  exists,  is  that  the  determinant  of  the 
coefficient  matrix  vanishes..    The  resulting  characteristic  equation  is 
called  the  dispersion  equation,  as  its  eigenvalues  determine  the  indices 
of  refraction  and  coefficients  of  attenuation.     That  is,  the  dispersion 
equation  is , 


€'        -n   2  -  n  2 
X  X       y  z 

€'        +  n     n 

y  x      y    x 

c'       +  n      n 
z  X         z      x 


c       +   n    n 

xy       x  y 


€        +  n    n 
Xz         X    z 


/  ?  2  e  '         n     n 

€        -n^-n  yz+yz 

,yy    k   *        Z        ,  2      2 

€         +nn  c      -  n   —  n 

zy  7     y  zz         y       x 


=    0   . 


(139) 
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4.21    Special  Cases 

As  an  indication  of  the  application  of  the  dispersion  equation,  some 

special  cases  will  be  considered.     For  this  purpose,  re-write  equations 

(68)  and  (103),  respectively,  in  the  following  forms.     For  E    , 

0 


-  y  x  H  =  j  toe 


1_  x  r  -j  x  y 

p2  _  y^        pTyl 

j     x  y        l-    x  r 


0 


0 

0 

1-     X 

r 


and  for  E 


-  y  x  H    =  j  0)  €„ 


Let 


y 


(140) 


0    -j     X  Y 

r2  -y2 

o  -l-  x_xl 

r2  -y2 
o         o 


0 

0 

1-     X 


(141) 


y  -  ya    =  j  k0n  =  j  k0n  (a      cos  0+  a    sin  Q  cos  <p+  a     sin  8  cos  cp) 


(142) 


Using  the  dielectric  tensor  from  equation  (140)  along  with  the  dispersion 
equation  (139),  and  also  substituting  from  equation  (142),  for  <b=  tt/2, 


1-  n     sin     9  -    TX  -  j  X  Y 

■p<^  _  y2   "p^  —  x2 


n     sin  0  cos  9 


j  X  Y 
r2  -  Y2 


i  -  xt  -  .rx    . 

f2    -  Y2 


n     sin  8 cos  8 


0 


and  for  <p  =  0, 


•  j  X  Y 


1_    rx       _  n2  sin  2  $ 

r2  -  y^  r 2  -  x2 

j  X  Y  2 

z^2      \r2    +  n     sin  8  cos  8  1-    TX 

r     "  Y  F2  -  Y2 


0 


1-n    cos     8  -    X 


+  n     sin  8  cos  8      0 


=  0     (143) 


-  n2 


nz  cos  L  0 


1-  X 


-n 


=0 


(144) 
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Equation  (143)  is  the  dispersion  equation  for  a  vertically  polarized 
wave  propagating  in  a  longitudinal  direction,  whereas  equation  (144) 
is  for  a  vertically  polarized  wave  propagating  transverse  to  the  geodesic 
field. 

Now,  in  equation  (143),  let  B  =  tr/2, 

-)XY  0 


1-  n    -    r  x 

r2  -  y2    r2  - x2 


j  x  Y 


^2         2 

r   -  x 


1  -  n    -  r  x 


2  _v2 


x< 


i-rx 


r2  -  x2 


Thus, 


(l-XJL 


and 


Y      _  y 


n    =  1  -      X 


2  2  2 

-nV  ~  (XY)       =      0, 


r±  y 


=  o 


(145) 


(146) 


(147) 


For  equation  (144)  with  6  =  tr/2, 


1  -  n      -   r  X 
IX  Y 


r2-Y2 


-j  X  Y 

r2-Y2 
i  -  rx 


r2  -y2 


1-       x 


-  n 


(148) 


one  solution  is , 


n    =  1  -    X 


(149) 


For  the  other  solution, 


(i  -    r  x       -n  )  a-  rx      )  -  (  x  y     )     =  o, 

tT2         2~ 

r    -  y 


r2-Y2 


2         2 

r    -yz 


n   = i -rx 


2       2 
X      Y 


(r2  -rx- y2)  (r2  -y2)    , 


or, 


n    =  1  - 


x_ 

r-    y 


r-x 
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(15  0) 


(151) 
(152) 


Solution  (147)  is  for  a  vertically  polarized  wave  in  the  longitudinal 
direction  whose  exciting  field  is  normal  to  the  geodesic  field,  and 
solutions  (149)  and  (152)  are  analogous  solutions  for  the  transverse  case. 

For  the  azimuth,  or  horizontal  component  of  the  exciting  field,  the 
dispersion  equation  f or  <p  =  77  /2 ,  that  is,  for  longitudinal  propagation, 
becomes , 


9  2 

-nz  sin     0      -j    XY 


2 

n     sin  9    cos  8 


r2-Y2 

l-n2  -     r  X 


r2-x2 


n     sin  9  cos  9 


l-n2  cos2  9  - X 


=    0 


(153) 


Expanding, 


j  X  Y  (0)    +  (1-  r    X 


2  2 

r    -  yz 


r2-x2 


-  n 


2) 


from  which  a  solution  is  obtained, 
n2  =  1  -    T     X 


-  sin  9  n     sin  9  cos  9 
cos  9    l-n2  cos2  9-     X 


=  0,     (154) 


Jl2     _    y2 


(155) 


For  the  other  factor, 


-(1-   X    )  sin  9  +  n     sin  9  cos     9-  n     sin  9  cos     9=  -(1-    X       sin  Q) , 

r  r 

which  is  constant.    Therefore,  solution  (155)  is  the  sole  solution  for  this 
case. 

For  transverse  propagation,  <p=  0,  and  the  exciting  field  is  parallel 
with  the  geodesic  field.     Hence,  in  this  case,  the  ordinary  solution 


holds,  namely 
n2  =  1-     X 


(156) 


5.    APPLETON  EQUATION  FOR  ARBITRARILY  ORIENTED  GEODESIC  FIELD 
5  .1    The  Appleton  Euqation 

The  Appleton  equation,  sometimes  referred  to  as  the  Appleton-Hartree 
equation,  is  an  equation  for  determining  the  complex  index  of  refraction  in 
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a  homogeneous,  anisotropic  ionosphere.     It  is  customarily  derived  by 
choosing  one  of  two  coordinate  planes  as  being  determined  by  the  y 
and  B0  vectors.     It  will  be  derived  herein  with  an  arbitrarily  oriented 
geodesic  field  so  that  a  greater  leeway  in  the  sfelection  of  the 
coordinate  axes  is  permissible  for  applications.    Refer  to  figure  8. 


*  #>> 


*£. 


Figure  8 
The  following  symbols  will  be  used, 

9  9  i 

W.   =  e  Bw    ,  jto      =  Ne       ,  T  =  1  -  j  v/w  ,  X  =  a).  z 

b  p 


m 


m  e„ 


/ 


CO 


f 


to, 
Y  =  _b_  a 

to 


2  2  2 

Y      +  Y        +  Y 

x       y         z 


(157) 


Y     =  Y  cos  6  ,  Y    =  Y  sin  9  cos  <p ,  Y    =  Y  sin  9  sin  <p  = 

x  y  z 


Y  cos  i/j  =  Y      , 


Ym  =  Y  sin  ip 


=  y/i- 


2  «     •     2 

sin     9  sin      <p  ,  e  B0     =  -  jy    . 


j  to  m 


The  geodesic  field  is  taken  as, 


BQ  =  B0  (a    cos  9+  a      sin  8  cos  <p  +  a       sin  9  sin  cp), 

x  y  z 


and  y  as, 


y     =az  y  =  j  k0    n 


(158) 
(159) 

(160) 
(161) 


5  .2    Lorentz  Conductivity  Tensor  for  An  Exciting  Wave 
The  exciting  field  will  be  postulated  to  vary  as  1       e 


jtot  -  y  .  r  . 


Hence,  Maxwell's  equations  may  be  written  in  the  form, 

y    .  H=0,y    .   c  0C  ' .  E  =  0,y  xE    =  j  w/x  0H  ,  -yxH=L  + 
j  toe0  E,  i  -  *  NeV  . 


(162) 
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Solving  for  i  , 

L  =  -  y  .    x  H    -  jo5€0  E  =  -  [j  05€0  E    +  _J yx(yxE)   ]  =  -  jwc  0/ 

(163) 


[  E  +   -  k  n      n  x  (ii  x  E  )   ]  . 


-  k 


Now, 


or, 


nx(nxE)  =  nxfhx(a      E     +a     E    +  a     E)]  =  n    a     x(a     E    - 
L  XXyyzz  zvyX 

a    E  )  , 

x  y 

nx(nxE)  =  -n     (a     E    +a     E      +  a     0). 

xx      y   y       z 


(164)' 
(165) 


Therefore , 

L  =  -j  wc0[E+nx(nx  E)J=  -  j  0)f  0  [  E  -  n     (a      E     +  a     E    +  a     0)], 

j  x    x      y   y     z 

(166) 


and  hence, 


n     -  1  0  0 

0  n2  -  1         0 

0  0-1 


(167) 


The  Lorentz  conductivity,  a  ,  is  defined  by, 
i  =  a   .  E 


(168) 


2  2         2 

Hence,  if  M     is  defined  by  M    =  n     -  1,  the  Lorentz  conductivity  tensor 


becomes , 
o  =  j  toe0 


n     -1        Q  O 

O        n     -1       O 

O  O  -1, 


=  J  <*>« 


M 
O 

O 


-• 

o 

\A2 

o 
o 

O 

-1 

1   U€~ 

M 


1    o 

o 

0     1 

o 

o   o 

-1 

M2J 

(169) 


A  corresponding  Lorentz  resiltivity  tensor,  p   ,  may  be  defined  by, 

p  .  I  =  p  .  a  .     E=F.E    =  E ,  (170) 

witirjl  being  the  identity  tensor. 
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The  resistivity  tensor  may  be  found  from  equation  (36),  which  will  be 
repeated  here  for  convenience, 


m     d  V      +  v  m    V+eV    xB0=    -eE. 
d  t 

In  the  steady  state  phasor  form,  equation  (36)  is, 

-eE    =  (  j  o)m  +  vm)    V    +  e  V  x  B0    0 

Using  equation  (160)  for  B0in  V  x  B0, 


VxB0     = 


or, 


a  a  a 

X  y  z 

V  V  V 

x  y  z 

cos  ©      sin  0cos  <p  sin  0sin<jp 


VxB  =  Bn[a    (V    sin0  sin  (p-V    sin0cos<p)  +  a     (V    cos  6  -  V 
9      X      y  z  y       z  j 

sin  0  sin  cp)  +  a     (V    sin  0  cos  (p-  V    cos  8  )  ]  . 

z      y  y 


(36) 


(171) 


(172) 


(173) 


Substituting  equation  (173)  into  equation  (171),  factoring    j  oj  m    from 
the  first  two  terms  of  (171)  and  using  the  symbol  T ,  the  resulting  equation 
is, 

-e  E  =  j  com    T(a     V+a     V.+a     V  )  +  e  Bn[a    (V     sin  0sin  <p- 

xxyy      zz  xy 

V    sin0cos<p)  +  a     (V    cos  0  -  V     sin0sin(p)  +  a    (V     sin0cos<n 
z  ^7yz  x  ZX 

-  V     cos  0  )  ]  . 

y 

» 

Substituting  symbols  (157)  into  the  re-arrangement  of  equation  (174)  as 
a  tensor  equation,  the  equation  may  be  written, 


E 
X 

r 

"iYz 

jV 

V 
X 

E 

y 

=    -  i  03  m 
e 

^Yz 

r 

-iY. 

• 

V 

y 

=  -1 
j  cj  e0x 

E 

-jY 

y 

iYx 

r 

V 
z 

r 

-jY            jY 
z           J    y 

— 

f 

X 

Ly 

jYz 

r      -rrx 

• 

/ 

"jYy 

jy        r 

X 

Lz 

— 

(175) 


with  the  velocity  matrix  converted  to  the  current  density  matrix  by 
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multiplying  and  dividing  by  -  Ne. 

From  equations  (170)  and  (175),  the  resistivity  tensor  may  be  taken 
to  be, 


=  -     1 


J    W€0  X 


r     -jy 


jy     r 

z 


JY 

y 
-JY 


-JY      jY 

y      x 


(176) 


Consider  equation  (17  0),  and  into  it,  substitute  from  equation  (188). 

The  result  is , 

E  =  p.T=p.a.    E  ,  T  .  E  -  p    .  <j     .  E    =  0 ,  (177) 

or, 

[n.g-l].E=0      ..  (178) 

Equations  (178)  constitute  a  set  of  the  linearly  homogeneous  equations 

for  the  three  components  of  E.     Hence,  for  a  non-trivial  solution  to  exist, 

p  .    a   -1       =0     .  (179) 

Multiplying  (169)  by  (176)  and  substituting  into  equation  (179),  the 

characteristic  equation  becomes, 
,6 


-    M 


X3 


r+  x 

M2 

")Yz 

]Y 

y 

M2 

iYz 

r  +  x 

"m2" 

"iY,* 

M2 

"jYy 

iYx 

r-x 

M2 

=  0 


(180) 


Expanding  the  determinant, 

(r+     X    J2     (T-X)    +jY    Y      Y-jYYY-Y2(r+      X    )  - 


M' 


x  y    z       x  y   z     y 


M< 


y2  (r+ jl.  )  -y2   (r  -x)  =  o, 

X  M2 


or, 


(r+    x    2  -  y  2 

M2* 


r  -x 


(  r  +  x__)  -  yt     =  o, 

M  L 


(181) 
(182) 


2  2  2 

in  which  substitutions  YT  =  Y    and  Ym    =  Y, .    +  Y        were  made. 

l      z         t        x       y 
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Solving  for  T+    X , 

M 

r  +  x    =1/2 

M2 

or,  transposing, 

YT 

r  -x 

_±/vT4 
1  (r-x2) 

+    4YL2 

X 

Y  2 

r-     t 
2  (r-  x) 

/     Y    4 

2 
+    4  YT 

M2 

1 4  (r-xj^ 

L 

(18  3) 


(184) 


Taking  the  reciprocal,  multiplying  by  X,  and  eliminating  M 

X 


™2         2     . 

M     =  n     -1  =  - 


r-  yt 


-I V    Y_  +Y 

2(F"X)         1  4-0^,2 


(185) 


Solving  for  n    and  substituting  from  equations  (157)  and  (159), 


n     =  1-        = 


o>p2  /  o>2 


2  2 

.    .     /        (aiJ^)    sin      0 
1-j  v/ui  -__b__; 


00 


oo 


b        cos      ib 
TT 
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2(1-0)  2/  toz-  jlTw) 


60,4/to     sin  j/) 


4(1-   60    Z7w2   -  JlVoo)2 
P 


(186) 


Equation  (186)  is  the  well  known  Appleton  equation  with  0  being  the  angle 
between  the  geodesic  field  and  the  direction  of  propagation.     It  is  not 
restricted  to  any  coordinate  system.    That  is ,  it  is  a  mathematical 
invariant. 

6 .     FARADAY  ROTATION 

6.1   Waves  Through  the  Ionosphere 

Waves  passing  through  the  ionosphere  are,  generally  speaking,  broken 
into  two  or  three  distinct  waves  having  distinct  indices  of  refraction,  and 
they  travel  by  distinct  ray  paths.     If  the  waves  are  attenuated,  they  are 
unlikely  to  recombine  into  linearly  polarized  waves. 
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However,  if  attenuation  is  negligible,  the  circularly  polarized 
waves  may  recombine  into  linearly  polarized  waves,  but  they  will  have 
their  plane  of  polarization  rotated  from  the  original  plane.    This  is 
referred  to  as  Faraday  rotation. 

To  examine  this  rotation  analytically,  consider  two  such  waves 
emerging  from  the  ionosphere  after  undergoing  different  path  length 
shifts  in  phases,  0  ,  0„   ,  respectively.    Write  the  phasor  equations, 

E  =  E0  e  "jk°r°  fe  ~jk°  Jn  ,  •  ds  ,    +  e    -jk0  J n2  .  ds"2  1 

Since  the  ionosphere  is  homogeneous  by  hypothesis, 


(187) 


E  =  E0  e  -jk°r°  (e  "jk°  nl  Sl    +  e  "jk°    n2  S2  )  (188) 


°r/    E  =  E0e-jk(r°    +     Vl+n2S2) 


3qJ^(nlSl-n2  s2)    + 


2  L      2 

eii^(n1s1-n2s2)   ] 

in  which  s]  and  s     are  the  respective  path  lengths. 
Upon  multiplying  and  dividing  by  2 , 


(189) 


E  =  2  E0  cos  [  ^     (nSi  -  n   "   J]  e'^^+  l/2  (Vl  +  n2S2)L         (190) 


2  112     2 


Thus,  the  resultant  of  the  shifted  vector  is  2E0cos   [k0  (n  s    -  n   s    ]  , 

and  k'0    (n^s    -  n2s2)  is  the  angle  of  the  resultant.     This  may  be  verified 
by  referring  to  figure  9 . 


<Pl  <f>2 

R  =  2  cos 


? 

<Pi  =     k0    n.s 
T  l 

^2=]ll   n2    S2 
2 

Figure  9 

From  equation  (190),  it  may  be  seen  that  the  phase  shift  of  the  field 

through  the  medium  is  determined  by  the  average  path  length  of  the  two 

ray  paths. 
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6 . 2    Polarization 


In  order  to  determine  a  measure  of  the  complex  polarization  of  a 
wave  propagating  within  a  homogeneous  anisotropic  ionosphere,  it  is 
desirable  to  find  the  ratios  of  the  electric  field  components  as  deter- 
mined by  equation  (179).    The  matrix  fcrm  is, 


M 


X 


r  +    x 

-  JY 

-JY 

.,2 

z 

V 

M 

M2 

iYz 

r+    x 

-n — 

+  JY 
X 

M2 

"iYy 

iYx 

r 


X 


=  0. 


L 

From  the  theory  of  linear  homogeneous  equations,  the  ratio  E    : 

A 

E     :  E     may  be  found  from  either  pair  of  the  three  equations,  given  by 
equation  (191),  by  omitting  the  first,  second,  and  third  columns  of 
coefficients,  respectively.    The  signs  are  alternately  plus,  minus, 
plus.    Thus, 


E    :  E 

x     y 


JY 


E     = 
z 


r  +  x 


M2 


JY. 


JY 

-*■ 

M2 

-(r-x) 


W 


JY 


-JY. 


JY 


# 


M 

-(r-x) 


w 


z 
-JY 


y 


r+  x_, 

JY. 


The  polarization  R  is  usually  defined  in  terms  of  the  ratio  of  two 
components  of  the  electric  field  normal  to  the  direction  of  ;phase  pro- 
pagation.   To  facilitate  the  algebra,  let, 

Y      =  0  ,  Y    =  YT    ,  and  Y    =  Y^ 
X  z         L  y        T 

This,  in  effect,  rotates  the  coordinate  axes  (figure  8)  such  that  B0 
lies  in  the  Yz  -  plane.     Now  write  the  ratio, 


R  ■ 


r+ x_ 

M2 

JY 

X 

jYx 

-or  ■ 

-x) 

jYz 

jYx 

-JY 

-(r- 

x) 

(.r-x)  (r+  x    )- y 

J  (r  -  X)  Y    +  Y    Y 

z      x    y 
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(191) 


(192) 


(193) 


(194) 


which,  upon  substituting  from  equations  (183)  and  (193),  becomes 


R  =  ^L 


2YL 


T  WYm  AJ2 


J 


±£ +  4Y  =    -J_0T+    X_  )        .       (195) 

(r-xr  ^         jyl        M7^ 


The  other  ratios  may  be  expressed, 


E-  2  Ez  2 

jYT    M  .     T~     =    °JYT    M       R.  (196) 


z       _ 


E 


y 


x 
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7.     CONCLUSION 

7.1    Equivalence  of  Two  Points  of  View 

For  the  purpose  of  yielding  a  better  insight  into  the  mechanics  involved 
when  a  wave  is  propagated  into  the  ionosphere,  the  principles  of  electron 
ballistics  were  applied  to  the  free  electrons.     In  particular,  equation  (36) 
was  solved  by  an  iterative  procedure.in  the  procedure .  for  deriving 
equation  (68),  the  latter  equation  yielding  the  dielectric  tensor. 

The  components  of  the  dielectric  tensor  were  used  in  equation  (139) 
for  determining  the  indices  of  refraction.    This  dispersion  equation  can  be 
formally  written, 

[n2  I  -  (n     +  V)]    .  E    =  0,  (197) 

in  which  an  index  tensor  n    is  introduced,  with  the  components  of  n 
being  defined  as , 

n  =  [n      ]=  [nj  n    ]    ,  i,  j  =  \>  Y,  z.  (198) 

Equation  (36)  was  also  used  as  a  key  equation  of  constraint  in 
deriving  the  Appleton  equation  (186).     Thus  equation  (36)  may  be  thought 
of  as  a  sort  of  common  denominator  between  the  two  procedures  for 
finding  the  complex  indices  of  refraction. 

In  fact,  equation  (68)  can  be  derived  much  more  compactly  by 
.'formalized  procedures.     For  this  purpose,  consider  the  corresponding 
complex  form  of  equation  (36), 

-eE=(ju)m+ym)V    +eVxB0  (199) 
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Referring  to  figure  5,  and  substituting 

-]Y  =  e  BB 
j  a>m 
equation  (199)  becomes, 

-e  E  =  j  a)  m  [(  1  =  j    v/u>)  V  -  j  Y   V  X  a     ]  =  j  w  m    x 


[r  V    +JYV    a 

x    y 


j  Y  V    a     ] 

y    x 


(200) 


(2  01) 


Hence, 
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o        or 
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I 
X 

I 
z 
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jw  e0X 


(202) 


From  equation  (202),  the  Lorentz  resistivity  tensor  n  may  be  written 
namely, 

r     -jy      o 

jy     r        o 

o 


=  -  i 


(203) 


Now,  the  complex  dielectric  tensor  is  given  by, 
=  '=   I     +  j  o>e0  (p)-l 


(2  04) 


Thus,  the  inverse  resistivity  tensor  may  be  found  by  customary  matrix 


algebra , 

p-1    =      -J0^    X 

r(rz  -x2 ) 


or,     =    , 
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r  x 
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(206) 
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Substituting  into  equation  (2  04)  and  introducing  symbols  (63)  and 


(64), 


l-    r  x 
x2  -  rz 

jXY 

x2  -  r2 

o  o 


■jXY 


X2  - 

l-  r  x 


x2  -  r 

o 


o 
o 

1-    x 


rj 


X 

-JY 

JY 

X 

0 

0 

(207) 


Substituting  c'     into 
curl  H    =  j  W€0  e    'E    , 
yields  precisely  equation  (68),  that  is,    _1 curl  H    = 

a     (E     X-j"E    Y)    +a     (jE    Y+E    X)    +  a     Z      . 

x    x  y  y       x        y  z 


(208) 


(68) 


The  above  derivation  tends  to  place  more  confidence  in  the  previous 
physical  interpretations.  It  also  serves  to  tie  the  iterative  procedure  to 
the  Appleton  equation. 

7  .2    Inhomoqeneous  Anisotropic  Ionosphere 

There  is  no  exact  mathematical  model  for  the  inhomogeneous  aniso- 
tropic ionosphere.     Many  statistical  measurements  and  mathematical 
interpretations  have  been  and  are  being  made.    A  vast  amount  of  lit- 
erature exists,  but  various  precise  studies  remain  to  be  made. 

A  detailed  study  of  the  inhomogeneous  ionosphere  was  entirely 
beyond  the  scope  of  the  time  and  facilities  available  for  the  pre- 
paration of  this  report. 
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